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Abstract 

We calculate the low temperature longitudinal ultrasonic attenuation rate as in clean d-wave 
superconductors. We consider the contribution of previously ignored processes involving the exci- 
tation of a pair of quasi-holes or quasi-particles. These processes, which are forbidden by energy 
conservation in conventional s-wave superconductors, have a finite phase space in d-wave super- 
conductors due to the presence of nodes in the gap which give rise to soft low-energy electronic 
excitations. We find the contribution to as from these processes to be proportional to T in the 
regime ksT <C Qua *C Ao, (ultra-low temperature regime) and to be proportional to 1/T in the 
region Qvp <C ksT <C Ao, (low temperature regime) where Q is the ultrasound wave- vector and 
Ao is the maximum gap amplitude. We explicitly evaluate these terms, for parameters appropri- 
ate to the cuprates, for Q along the nodal and the antinodal directions and compare it with the 
contribution from processes considered earlier(I.Vekhter et al Phys. i?eu.B59, 7123(1999)). In the 
ultra-low temperature regime, the processes considered by us make a contribution which is smaller 
by about a factor of 10 for Q along the nodal direction, while along the antinodal direction it is 
larger by a factor of 100 or so. In the low temperature regime on the other hand the contribu- 
tion made by these terms is small. However taken together with the original terms we describe a 
possible way to evaluate the parameter vf/va.- 

PACS numbers: 74.20.Rp; 74.25.Ld 



1 



Ultrasonic attenuation has been an important probe in the understanding of BCS su- 
perconductors. One of the earliest verifications of the temperature dependence of the gap§, 
predicted by BCS theory, was made by measuring the ultrasonic attenuation co-efficient 
(UAC). Thus, it is expected to shed light on the nature of the superconducting state in the 
cuprates as well. The calculation of the low temperature attenuation rate in d-wave super- 
conductors has been considered earlier, both within the clean limiti'iS and in the opposite 
hydrodynamic limit@i. Even earlier ultrasonic attenuation in anisotropic superconductors 
has received attention both in the cleanB and in the dirtyi'i'00 limits. On the experimental 
front, there exists just the one experiment^ but the results of our work, which have been 
obtained in the limit Ql ^> 1, that is the clean limit (where Q is the ultrasound wave- vector 
and I is the electronic mean free path) are inapplicable there. Assuming the ultrasound fre- 
quency to be 150 MHz, which lies in the range in which typical ultrasound experiments@0 
are done, and the sound velocity^ to be 4 x 10 5 m/s, the clean limitH corresponds to the 
condition I ^> 4/xm. Thus our calculations are applicable only to the cleanest sample 
High purity YBCO crystals have been fabricated in BaZrO^ cruciblesElhaving / ~ Afim, 
giving rise to hope that the clean limit may be experimentally realizable. 

There exists evidence^ that although the normal state in the cuprates remains a puzzle 
in the absence of well-defined quasi particles, the superconducting state is well-described by 
BCS-like quasi particles suitably adapted to the symmetry of the order parameter in the 
system. Further, experimental evidence largely points towards a d-wave order parameter^! 
in the high T c materials. The layered nature of the cuprate superconductors together with 
strong electronic correlation leads to a highly anisotropic electronic transport in the normal 
state, even though the superconducting state has three-dimensional phase coherence due to 
inter planar Josephson coupling. However most superconducting properties can be described 
by assuming 2-d BCS quasi-particles that have an energy gap with a d x 2_ y 2 symmetry. For 
the purposes of this paper we assume that the d-wave superconducting quasi-particle energy 
dispersion to be phenomenological. 

The quasi particles have a gap function that has nodes on the Fermi surface at the 
points k x = iky. At low energies, the relevant low energy electronic excitations, are the 
Dirac quasiparicles near these nodes. These nodal quasi-particles have energies that are 



well described by = yvpkf + v\k\. This Dirac spectrum is easily obtained^ by a linear 
expansion of the band energy and gap about the nodes of the single particle excitation 
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spectrum. Here k\ and ki are the momentum components normal and tangential to the Fermi 
surface at the nodal points. While it is possible to obtain such a quasiparticle spectrum 
from a suitably chosen microscopic HamiltonianIS in this work we treat Vp and v& as 
phenomenological parameters whose values should be read from experiments. Because of 
the dominance of the band energy scale W > A , where W is the band width and A is 
the maximum amplitude of the superconducting gap, ^ is large and is approximately 10 
for parameters relevant to the cuprates. 

The UAC is the inverse of the length obtained by a product of the sound velocity and the 
phonon life-time and sets the scale over which an ultrasonic wave decays while propagating 
in the solid. The temperature dependent variation of the sound velocity is very weak in 
the temperature range of interest to us. We therefore calculate only the dependence of the 
phonon life-time on temperature. In general the damping of sound waves in the lattice can 
arise due to various excitations in the system. Here we focus exclusively on the damping of 
phonons due to their interaction with electrons. We ignore the contribution of anharmonic 
terms in the phonon Hamiltonian, as these are unaffected by the onset of superconductiv- 
ity. Thus, assuming the usual electron-phonon coupling, the UAC is proportional to the 
imaginary part of the electronic density- density correlation function. 

A calculation of the low temperature UAC in d-wave superconductors was first performed 
by Vekhter et. al.i and generalized to the orthorhombic situation by Wu and CarbottJ. 
They obtained an attenuation rate which is proportional to the temperature at T Ao. To 
obtain a result that is independent of the poorly known electron-phonon coupling constant, 
they normalise their result by dividing the low temperature attenuation rate (as(T)) by the 
attenuation rate at T c . In order to obtain a result analytically for as{T c ), it is necessary to 
assume a parabolic band dispersion. It is by now well-established0 that the Fermi surface 
in the cuprates deviates strongly from circular symmetry. Thus there is an unknown , Fermi 
surface sensitive wave vector dependence in as{T c ), which has been missed. A fundamental 
drawback of their calculation is that the use of a quasiparticle description at T c is itself 
questionable given the anomalous nature of the normal state in the cuprates. 

Apart from these questions, all previous studies of the UAC in references Refs. (2,3,4,7) 
ignore the contribution of processes corresponding to the creation of a pair of quasi particles 
or a pair of quasi holes. In conventional s-wave superconductors, where the condition Q <C A 
is satisfied, energy conservation forbids these processes. However, in d-wave superconductors 
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the nodes in the energy gap imply that there is a finite phase space for these processes. In 
Ref.(3) it is claimed that phase space restrictions make the contributions accruing to the 
UAC from here negligibly small. In this paper we explicitly evaluate the contributions of 
these processes in the low temperature regime. At ultra low temperatures (hsT <C Qv& <C 
A ) where Q is the ultrasound wave-vector, vp is the Fermi velocity at the nodes and A 
is the maximum amplitude of the superconducting gap, we find that there is a term in as 
linear in T from these processes in addition to a similar contribution evaluated earlier by 
Vekhter et. al. However, the co-efficient of the linear in T term, in this case is found to be, 
for parameters appropriate to the cuprates, larger by a factor of 100 for sound wave along 
the antinodal direction. In the nodal direction these new terms provide corrections that are 
smaller by a factor of 10. 

The other regime we consider is Qvf <C ksT <C A and here we find the contribution 
to the UAC from these previously neglected processes is proportional to 1/T. Using the 
ratios of the co-efficients of the 1/T contribution calculated by us and the linear in T 
contribution evaluated in Ref.(2) for two different ultrasound wave- vectors Qi and Q 2 it is 
possible to determine the ratio vf/va which is a basic parameter in the description of the 
superconducting state in the cuprates. 

The imaginary part of the density-density response function is given by x"(Q,f2) = 
X a (Q, n) + xl(Q, n). Here X " a (Q, V) = XI - X2 and Xb (Q, = X3 - X4, where 

XI = ^ k [n(E k ) - n(£ k+Q )]« Q 4 " 4 ^^ W" + E * ~ ^+Q)> C 1 ) 

X2 = ^ k [n(E k ) - n(E k+Q )](vl +Q vl - ^^ )S(Q -E k + £ k+Q ) (2) 

*3 = ^E k [l - n(E k ) - n(E k+Q )](ul +Q vl + ^£^ )W E k+Q ), (3) 

and 

= ^S k [l - n(E k ) - n(E k+Q )](vl +Q ul + ^g| )a(n + E k + E k+Q ). (4) 

To proceed further, we retain terms only to the leading order in Q and f2. To this order 
x" a reduces to 

£ = ^Z k n>(E k )5(Q.v g )^ (5) 



where v g = V k -E k . 
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This term has been calculated earlier by the authors of Ref.(2). For k B T <C A , it yields 

-Aln2Vta 2 k B T f3 2 A 2 
Xa ~ 2^ ^^}{a 2 + /? 2 ) 3 / 2 + (A 2 + r] 2 )^ ^ 

Here a = vpQ cos(7r/4 — 9), (3 = t>A<3 sin(7r/4 — 9), rj = vpQ sin(7r/4 — 9) and A = 
t>AQcos(7r/4 — 9). We now turn our attention to the remaining contribution to x" which 
comes from the terms in Eqs. (3) and (4). In conventional s-wave superconductors we have 
Q <C 2A and so the energy conserving delta function in Eqs. (3) and (4) can never be satis- 
fied. Thus xKQ)^) — in that case. However, as mentioned earlier, the presence of nodes 
in the gap function leads to a finite contribution from x&(Q? ^) m d-wave superconductors. 
It has always been assumed in previous work@Ji0Bil'0 that the phase space available for the 
processes described by this term (creation or annihiliation of a pair of quasi particles) would 
be small and so its contribution has been ignored. In this paper we explicitly evaluate this 
term. We find that its contribution is large and cannot be ignored at ultra low(T <C lmK) 
temperatures. In the low temperature regime (lmK <T< A //cb) although the relative 
size of the terms evaluated by us is small its different (1/T) temperature dependence enables 
the determination of the parameter vp/v& by measuring the attenuation rate at 2 different 
angles. These are the principal results of this paper. 

We now present our calculation of x&(Q? AVe work in units of lattice constant a = 1. 
To linear order in Q, we find 

x b (Q, n) = ^x k n'(£ k )K +Q V 2 + WV + ft kZ i k+Q )5(E k + £ k+Q ) (7) 

To leading order in Q, Xb" further reduces to 

xJ(Q,«) = ^ w nn'(E k )5(2E k + Q.v g ) (^) 2 (8) 

Linearizing the band energies and the gap function at the nodes0 which lie at the in- 
tersection of the lines k x = ±k y with the Fermi surface, we find that £ k ~ vpki, A k = 



ua&2 and Ek = \jv 2 F k 2 + v\k\ where k\ and k 2 are the components of the momenta normal 
and tangential to the Fermi surface at the nodes. We take Q = Qcos9x + Qsin9y. Going 
over to polar coordinates (f^/ci = pcoscf) and v a&2 = psinip) we find that for the node near 
(7r/2, 7r/2), Vg.Q = acoscj) + fisincf). 

In comparison to the corresponding contribution for Xa" an extra term 2p appears in 
the argument of the delta function in Eqs. (3) and (4). Finally we obtain for the node near 



(n/2,n/2), 

„ Q /"Pc /"27T 

X 61 (Q,fi) = - / dppn'(p) / d(f)5(2p + acos(f) + (3sin(j))sin 2 4> (9) 

2~/'/- ''a •' u ./ii 

where p c = ^/twjPua is a cutoff chosen to preserve the volume of the Brillouin zone. We 
perform the integral first. The condition that the zeros of 2p + acos<p + (5sin(f) are real 
together with the requirement that cos 2 0q < 1 (where 0q are the 2 independent zeroes in 



the range < <p < 2n) leads to the constraint p < — — ^ — = p c . Assuming Q ps 400cm 
and va/vf ~ ^ ~ 1/10 for the cuprates we find that p c <C p c . Thus the contribution to 
Xb" from this node can be written as 

8p 2 (/3 2 -a 2 ) , 2a 2 

" (r\ n\ 11 f A '( \ (" 2 +/3 2 ) 2 a 2 +/3 2 / ln x 

X - (Q ' fi) = 2^io ^ (p) (^ + ^-4p 2 )V 2 ^ 
The change of the upper cutoff from p c to p c is the phase space reduction referred to by 

the authors of Ref.3. 

We evaluate the expression in Eq.(10) in two regimes. &)k B T <C Qva <C A and b)Qv F <C 
k B T <C A . The maximum value of p c is p c (max) = vpQ/2 and minimum value of p c is 
p c (min) = vaQ/2. We have taken the temperature regimes to be either a) much less than 
p c (min) or b) much greater than p c (max). 

We approximate the derivative of the Fermi distribution so as to simplify the integration 
and render it doable analytically. We notice that n'(p) ~ — l/Ak B T for p/2k B T <C 1 and 
n '(p) ~ ~ e Jy gT for p/2k B T 3> 1. We therefore make the interpolating approximation 
n'(p) — l/Ak B T for p < 2k B T and n'(p) ps ~ £ '' k "^ BT for p > 2k B T. Our approximation is 
asymptotically exact and is therefore expected to provide reasonable answers in the 2 limits 
{k B T <C p c <C A and p c <C k B T <C A ) considered by us. 

We first consider the ultra low temperature regime k B T <C Qva- For typical parameters 
appropriate to the cuprates and assuming Q ps 400cm~l this condition implies k B T <C lmK. 
In this case x'bi can be split into 2 pieces and can be written as 

x?i(Q,n) = (— i/4/c b t) — — — [i i + i 2 ] (ii) 

TTVfVa 



where 



and 



,2fc s T p ^2(^2 _ Q 2) ft 2 

1_ 7o + - V (a 2 + /3 2 ) 2 + a 2 + /3 2) ( } 

rPc pe^ 4p 2 Q3 2 -a 2 ) a 2 

2fe B T P y/d* + ^ " V (« 2 + /3 2 ) 2 « 2 + /3 2J 1 ' 
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The integral in l\ is elementary. To evaluate I2 we further make the approximation 
J^dxe~ x f(x) ~ e~ l J™ OT ('+ 1 >' U ) dxf(x) (for u > I > 0) and where f(x) is a slowly varying 
function compared to e~ x . We similarly evaluate the contributions of the other three nodes 
as well. The final answer is a long expression and in the interest of clarity and brevity 
we only present the leading term in an expansion in the small parameter , kBT . To this 
leading order we find 

„ —ksTVt _ 2 a 2 rj 2 

Xb = nv F v A (1 + 56 )[ {a 2 + (3 2 fl 2 + ( V 2 + A 2 ) 3 ^ ( > 

It is instructive to compare this result with the processes considered earlier in Eq.(6). We find 
that upto factors of 0(1), the result is the same with a <-> /3 and A <-> rj. This change comes 
from the different factor ^ ~ sin 2 (f) in Eq. (8) as compared to J^- ~ cos 2 (f) in Eq. (5). Also 

k k 

notice that at ultralow temperatures the dimensionless cutoff to the integral p c jksT ^> 1 so 
that the phase-space reduction due to the replacement of p c by p c is ineffective as most of the 
contribution to the integral comes from the region where p/ksT <C p c //c#T. To get an idea 
of the relative size of the contribution of the processes considered by us with respect to those 
considered earlier by Vekhter et. al.i we have evaluated % in the ultra low temperature 

Xa 

regime for Q directed along a node and an antinode respectively. In the nodal case we have 
9 = 7r/4, a = vpQ, (3 = rj = and A = v&Q. Then from Eqs. (6) and (14) we find that 

Xl _ (l + 5e- 2 ) v A 

X'L 2/n2 v F 1 ' 

Thus the contribution from these terms is smaller by a factor of v^/vp ~ 1/10. For the 
antinodal direction we have 9 = 0, a = 77 = vfQ/Vz, (3 = X = vaQ/\/2. Putting these 
values in Eqs. (6) and (14) we have 

Xl _ (1 + he- 2 ) v\ 

Xl 2/n2 v 2 A { ) 

Thus the ratio in this case is approximately 100 and so in this case it is the processes 
calculated by us that make the dominant contribution. 

We now proceed to the evaluation of x'l m the limit Qvp ksT C Ao. In this case, we 
have n'(p) ~ — l/AksT over the entire range of p integration from to p c . Thus evaluating 
and summing the contribution to x'l fr° m eacn of the nodes we arrive at the result 

» g 1 ( q2+2 P 2 x2+2r ? \ <m 

Xb 2A7Tv F v A k B T\a 2 + (3 2 y/ 2 {X 2 + r] 2 y/ 2) 1 ' 



In this case the relative size of x'l with respect to x" a * s veT Y small 0((QvF/kBT) 2 ). 
However we now show how our result may be used to extract the parameter The total 
imaginary part of the density-density correlation function can be written as 

x'LaiiQ, n) = xl(Q, n) + x?(Q, n) = 7i(Q, n)T + (is) 

where 71 and 72 are coefficients whose values are given in Eqs.(6) and (17). Thus the UAC 
is given by 

a s (Q, T) = M(Q)[ 7l (Q, Q)T + (19) 

where M(Q) is a factor dependent on the electron-phonon coupling constant, the sound 
velocity and the ultrasound wave-vector. A measurement of as(Q,T) and a comparision 
of the two terms contributing to it would enable an extraction of the parameter X(Q) = 
—(q 2 ^) - The poorly known factor M(Q) exactly cancels out. On taking the ratio of X(Q) 
evaluated at two different wave vectors we obtain an expression that depends only on the 
parameter ^ and the wave vectors Qi and Q2 respectively. For the sake of illustration we 
have computed X node {(^-i) / X antinode {Q l 2) for Q x along the nodal direction and Q 2 along the 
antinodal direction. We find 

X node {^)/X antmode ^ 2 ) - {{vF/vA)2 + mvF/vA y + 2) (20) 

Thus the measurement of as(Q,T) in the regime (10 _3 i^ <C ksT <^ Aq) provides an 
useful method to extract vf/va without use of as(Q,T c ) whose calculation necessitates 
assumptions about the existence of quasi-particles at T c with a circular Fermi surface. 

Finally let us recapitulate the main points of this paper. We have calculated the UAC 
for a d-wave superconductor. We have included the effects of processes involving creation 
or annihilation of a pair of quasi-particles which are non zero because of the soft excitation 
spectrum of a d-wave superconductor. In the ultra-low temperature regime we find that 
these processes make a substantial contribution, which is linear in T, to the UAC. In the 
low temperature regime we find a small 1/T contribution to the UAC. We indicate how 
this can be used to extract the parameter vf/va- A drawback of our work is that we 
have carried out calculations in the clean limit which makes it inapplicable to all but the 
cleanest of samples. Yet another drawback is that we assume the existence of quasi-particles 
of spectral weight unity, thus neglecting all effects coming from the incoherent part of the 
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spectral function seen in ARPES 
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